Orientation-dependent Casimir force arising from highly anisotropic crystals: 

application to Bi2Sr2CaCu208+5 
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We calculate the Casimir interaction between parallel planar crystals of Au and the anisotropic 
cuprate superconductor Bi2Sr2CaCu208+i (BSCCO), with BSCCO's optical axis either parallel or 
perpendicular to the crystal surface, using suitable generahzations of the Lifshitz theory. We find 
that the strong anisotropy of the BSCCO permittivity gives rise to a difference in the Casimir force 
between the two orientations of the optical axis, which depends on distance and is of order 10-20% 
at the experimentally accessible separations 10 to 5000 nm. 
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The Casimir force, an attraction between conductors 
that arises from the quantum mechanics of electromag- 
netism in the vacuum, has been of interest to funda- 
mental physics since its first description by Casimir 
Recently, however, interest has grown in "tailoring" the 
Casimir force by adjusting material properties like film 
reflectivity [1| , thickness [1] , and carrier density [1] . Some 
recent proposals suggest that it might even be possible 
to reverse the sign of the Casimir force between negative- 
index materials ^] or other materials with unusual elec- 
tromagnetic properties @, which could be realized in 
properly designed metamaterials [7|- Since most, if not 
all, designs for such metamaterials are anisotropic, it is 
necessary to consider what effect this anisotropy has on 
the Casimir interaction. 

Here we show that strong anisotropy in the dielec- 
tric permittivity alone has a significant influence on the 
Casimir interaction. We calculate the Casimir force be- 
tween crystals of Au and the cuprate superconductor 
Bi2Sr2CaCu208+5 (BSCCO); when the optical axis of 
the BSCCO crystal is switched between being parallel 
and perpendicular to the crystal surface, the Casimir 
force changes by a distance-dependent amount, up to 
25%, at the experimentally accessible separations 10 to 
5000 nm. BSCCO is well studied as a high-temperature 
cuprate superconductor, though superconductivity is not 
at issue here. Rather, BSCCO is an exemplary material 
for the present study because it shows great anisotropy 
over a very wide frequency range, yet is a homogeneous 
single crystal and can be modeled with its dielectric 
permittivity without having to consider its microscopic 
structure, as must be done at high frequencies with arti- 
ficial metamaterials. 

Some relevant theoretical work has been done already. 
In the van der Waals force was calculated between 
anisotropic plates in the non-retarded limit, which is only 
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applicable to separations of tens of A or less. In [9*], the 
Casimir-Polder force was calculated between an atom and 
a uniaxial crystal surface in one orientation (called the 
perpendicular cleave in the present work; see below) , and 
in [lO| the Casimir force between an isotropic surface and 
a uniaxial surface in the same orientation, both papers 
dealing primarily with graphite and graphene. Recently, 
[m studied specifically the Casimir force between a metal 
and an anisotropic metamaterial, though again only in 
one orientation (the perpendicular cleave), and at only 
one separation. In [12l . Il3l | , the Casimir-Lifshitz interac- 
tion between two uniaxial surfaces was analyzed in a dif- 
ferent orientation (called the parallel cleave in the present 
work) with the goal of calculating a quantum electrody- 
namic torque, but the relation of anisotropy to the dis- 
tance dependence of the force was largely unexplored. In 
the following, we show that in fact this distance depen- 
dence can be made to vary significantly by changing the 
orientation of the crystals. 

A uniaxial crystal is characterized by a dielectric per- 
mittivity that has the value e|| (w) for electric fields polar- 
ized along the so-called optical axis, and a different value 
for electric fields polarized perpendicular to the 
optical axis. The permittivity depends on the angular 
frequency uj of the field, but the direction of the optical 
axis is the same for all frequencies in a uniaxial crys- 
tal. Such crystals generally are birefringent, and waves 
with the same frequency in the same direction will have a 
wave number (or wavelength) that depends on their po- 
larization. Therefore if a cavity is composed of uniaxial 
crystal surfaces, the electromagnetic normal modes will 
have different dispersion relations for the two principal 
polarizations, which affects the zero-point energy of the 
modes, and in turn the Casimir interaction. 

In the present paper we calculate the Casimir force be- 
tween uniaxial crystal surfaces in two orientations (shown 
schematically in Fig. [1]). The first orientation, which we 
call the perpendicular cleave, has the optical axis per- 
pendicular to the surfaces. The second, which we call 
the parallel cleave, has the optical axis in the plane of 
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(i^) is the reflection ampUtude at an in- 



terface between vacuum and material j for a wave whose 
transverse component of the wave vector is k, imaginary 
frequency is i^, and polarization p £ {TE, TM}, and 
Kg = + ^ l(? ■ The reflection amplitudes are well 
defined for the perpendicular cleave because TE and TM 
are the principal polarizations for waves traveling in the 
uniaxial crystals (but not for other orientations of the 
optical axes). If the crystal- vacuum interfaces lie in the 
xy plane, the optical axis points in the z direction and 
the permittivity in material j is the tensor 



FIG. 1: The Casimir force between Au and BSCCO depends 
on the orientation of the BSCCO optical axis. Solid curve 
is fperp, the force per unit area for the perpendicular cleave, 
dashed curve is fpar for the parallel cleave, both normalized to 
the result for ideal metal plates fideai = —t^ hc/2'lOd^. The 
insets show the orientation of the conducting copper oxide 
planes in BSCCO. The optical axis is perpendicular to these 
planes. 



the surfaces. (We note that the parallel cleave can give 
rise to a QED torque as shown in (T^. ITsj.) 

We consider first the perpendicular cleave. (A similar 
derivation is used in Using Eqs. 86 and 90 from 

p^ . we calculate the Casimir force per unit area /(d) for 
parallel plates, made of media 1 and 2, separated by a 
distance d of vacuum: 



f{d) = 2hJ2 
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(We suppress the oj dependence of all permittivities for 
readability.) The TE and TM reflection amplitudes are 
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where kq = a/Zs^ — lo"^ /c^, kj = ^k"^ — e^^w^/c^, and 
Aj — -\/fc^ej_L/ej|| — €j\_LO^ jd^. Thus for this cleave 
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where we have made the substitution uj 

so Ko = \fW^\~^^ , Kj = ^/WTTjjJ^Jc? , Xj = 

^fc-^ej_L/ej|| 4- ejj_^2/c2, and the permittivities ej = 
are evaluated at imaginary frequency i^^ for ma- 
terials j = 1,2. For isotropic materials, e 
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and we recover the standard Lifshitz formula 



In the parallel cleave orientation, we cannot use the 
same method because the principal polarizations in the 
uniaxial crystals do not coincide with the TE and TM 
polarizations. Therefore, the normal mode frequencies 
must be calculated for the entire cavity at once, keeping 
track of all three components of both electric and mag- 
netic fields at both interfaces. This straightforward but 
quite cumbersome calculation is performed in (and 
the result is reported also in [l3|), yielding the free energy 
as a function of temperature and relative angle between 



the optical axes of the two materials. In the case at hand, 
we choose BSCCO as material 1, Au as material 2, and 
vacuum as the intervening medium 3, and also take the 
zero-temperature limit of the expression. (For the more 
general and more unwieldy form with anisotropic mate- 
rials 1 and 2, arbitrary intervening medium, and nonzero 
temperature, see Eqs. 2.1-2.11 of ^].) In this case the 
permittivity tensors are 




£2 = £2 




(5) 



With medium 2 and 3 isotropic, we obtain a greatly sim- 
plified expression for the Casimir potential energy per 
unit area Upar{d): 
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where 
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Ki = - (1 - ei||/eii)fc2 cos2 + eiii^Vc^, (9) 

KO = yp+ev^- 

As before, k = (fccos0, fcsini/)) is the transverse compo- 
nent of the wave vector with fc, <j) its polar coordinates, 
and the dielectric functions tj = ^jiiS,) are evaluated at 
imaginary frequency i^. The Casimir force per unit area 
is fpar{d) = —dupar/dd. If material 1 is isotropic, as well 
as material 2, ei± = ei|| and we recover the standard Lif- 
shitz formula. 

We now evaluate the Casimir force for our 
example case, taking the cuprate superconductor 
Bi2Sr2CaCu208+5 (BSCCO) for material 1 and Au for 
material 2. We assume the BSCCO to be optimally 
doped, and consider its normal-state permittivity. Like 
other cuprate superconductors, BSCCO is composed of 
layers of conductive copper oxide planes separated by in- 
sulating oxides; the conductivity is much higher in the 
copper oxide plane than perpendicular to it, and the 
permittivity is likewise anisotropic. (There is also a 
much smaller anisotropy between the two principal axes 
within the copper oxide planes p3|, which we neglect.) 
Due to the extreme anisotropy of BSCCO, the Casimir 
force differs in the two orientations by up to 10% for 
10 nm < d < 250 nm, and by up to 25% at larger sepa- 
rations, as shown in Fig. (TJ forces are normalized to the 
ideal metal result Q f^deal = -TT^nc/240d\ The ratio 
fpar/fperp IS shown in Fig. [21 For d > 25 nm, the force 
is stronger for the perpendicular cleave than the parallel 
cleave. For d < 25 nm, the reverse is true. 

To model the permittivity function of BSCCO, we use 



1.05 - 
1.00 - 
0.95 - 
0.90 - 
0.85 - 
0.80 - 
0.75 



r 

10-^ 



-1 1— 

4 



I I 1 1 1 1 1 — r 

6 8 2 4 

10-^ 

Separation [m] 



G 8 



10' 



FIG. 2: The ratio of Casimir forces fpar/fperp between 
BSCCO and Au, for the two orientations of the optical axis. 
The two forces are most different at large separations, where 
long wavelength permittivity dominates and the anisotropy 
of BSCCO is most pronounced. 



multi-oscillator models where 
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We model the permittivity of Au using a Drude compo- 
nent plus two interband transitions in the critical point 
transition model [l^ (since Au is not well described by 
a small number of Lorentz oscillators): 



2 

Cm(^m + 7m + C) ^^^-^ 



The parameters are shown in Table U (based on [T^ for 
Au, on [1^ for ej_ of BSCCO, and on our fits to reflectiv- 
ity and permittivity data [U for e|| of BSCCO). The 
anisotropy of BSCCO in our model is shown in Fig. O 
the top panel shows a much higher reflectivity (at nor- 
mal incidence) for waves polarized perpendicular to the 
optical axis, compared with waves polarized parallel to 
the optical axis, over a wide range of frequencies; the 
bottom panel shows the same information presented as 
permittivity versus imaginary frequency. 

To perform the computations, we use Mathematica 
6 to numerically integrate the expressions for fperp and 
fpar- To control a spurious contribution from roundoff 
error, we include an upper cutoff in the integration for 
the transverse wave vector and frequency of 10^ m~^ and 
c X 10^ m~^, respectively. Increasing the cutoffs by a 
factor of 10 changes our results insignificantly, by less 
than 1 part in 10^, at all separations. 

The difference in force between the two orientations fol- 
lows our expectations, at least for d > 25 nm, if we imag- 
ine BSCCO 's copper oxide planes as a stack of partially 
reflective mirrors. In the perpendicular cleave, these 
"mirrors" are parallel to the crystal surface and, with 
the Au surface, constitute a confining Fabry-Perot-like 
cavity for all modes. In the parallel cleave, the "mirrors" 
are perpendicular to the crystal surface and resemble a 
linear polarizer, and effectively confine modes of only one 
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FIG. 3: Optical properties of BSCCO. Waves with electric 
field perpendicular to the optical axis meet with a very con- 
ductive and reflective material, compared with waves polar- 
ized parallel to the optical axis. Top panel: reflectivity at 
normal incidence vs. real frequency for the two polarizations. 
Bottom panel: permittivity vs. imaginary frequency for the 
two polarizations. 
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TABLE I: Parameters used in models of the permittivities of 
BSCCO and Au. Resonance frequencies uurn, widths 7™, and 
strengths Am and Cm are in rad/s, and too is unitless. Powers 
of 10 are bracketed, e.g. 2[4] = 2 x 10*. See Egs. [TOlandfTT] 



polarization. We should expect the first case to have a 
stronger Casimir force than the second, and this indeed 
happens. At very small separations, the high frequency 
(visible and UV) optical properties become increasingly 
relevant, and at these frequencies the reflectivities for the 
different polarizations are similar. The analogy to mir- 
rors therefore breaks down at very small separations. 

The calculation of the Casimir force depends on the 
permittivity over all frequencies, and the incomplete 
knowledge of the permittivity function is a major source 



of uncertainty. Even sample-to-sample variations in the 
optical properties of nominally the same material can 
change the force by 5% or more [12] ■ This is the most 
significant uncertainty in our calculation for Au and 
BSCCO. Since the reflectivity of BSCCO over the re- 
ported frequency range, 30-30000 cm~^, is quite differ- 
ent from zero for both polarizations ^2Q], we can infer 
that there are resonances at higher frequency, but can- 
not know for certain their frequency or number based on 
the data published so far. 

For this reason, we add to £±{'-o) a single oscillator term 
A'jjp/{ujjjp ~ up' — iw^hf) with center frequency ijJhf 
much higher than the highest frequency for which we 
have data. This contributes a constant term Aj^p/uj'jjp 
for all uj <€. lohf- We chose uhf = 1-88 x 10^^ rad/s, 
corresponding to a wavelength of 100 nm. A lower ujhf 
would conflict with existing optical data. Increasing ujhf 
(while holding Aj^p/uif^p constant) does not have a sig- 
nificant effect on our conclusions: raising it by a factor 
of 10 changes the the ratio fperp/fpar by less than 1% at 
separations d < 100 nm (though the overall magntitude 
of fperp and fpar changcs by a few percent), and has even 
less effect at larger separations. Further increasing luhf 
does not further change the Casimir interaction. 

To choose our model for BSCCO's e||(w), we used the 
published reflectivity data [13] as well as un pub lished 
permittivity data covering 0.1 to 4000 cm~^ [2l[. We 
obtained a satisfactory fit to the optical data using a 
highest resonance at 33,000 cm~^. 

Regardless of these theoretical uncertainties, this cal- 
culation serves as a proof of principle that realistic ma- 
terial systems can show a significant change in Casimir 
force when their optical axis orientation is changed, even 
if the calculation is not quantitatively exact for the 
Au/BSCCO system. In any case, we expect our conclu- 
sions to hold qualitatively for this system, as they stem 
from the established anisotropic permittivity of BSCCO 
in the dc to infrared frequency range. 

In conclusion, we have shown that strong anisotropy 
in the permittivity of a uniaxial crystal can give rise to a 
Casimir force that depends significantly on the orienta- 
tion of the optical axis. We find that the Casimir force 
between Au and BSCCO plates varies by an amount 
of order 10-20 % when the optical axis of BSCCO lies 
perpendicular versus parallel to the crystal surface. We 
expect a similar effect for other materials with similar 
anisotropy. Such an effect represents a way to vary the 
strength of the Casimir force between two given mate- 
rials merely by changing their relative orientation. This 
adds another to the small number of "handles" on the 
Casimir force, which is usually determined only by ma- 
terial properties and the distance between the plates. 
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